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(i) 


ABSTRACT 


The object of this thesis is to determine and study the 
approximate gravitational field of an infinitismally thin rotating 


disk. 


The metric is reviewed and the linearized gravitational 
field obtained. A survey is made of hypersurfaces in general rela- 
tivity and the results are used to determine the energy tensor of 
the disk. An approach is outlined which enables one to determine 
the second order corrections to the field and energy tensor. The 
properties of the metric on the disk and in the far field limit are 
discussed. A short account of the gravitational red shift and the 


Lense-Thirring effect.is given. 
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CHAPTER I 


The Linearized Gravitational Field of a Rotating Disk 


in introduction. 


The purpose of this dissertation is to examine the approximate 
gravitational field of a thin rotating disk of gravitating particles. 
In this section we shall outline the arguments which lead to the most 
general metric which determines a stationary axially symmetric gravita- 


tional field. 


* 
In general the line element has the form 


du“¥e Baye dx” He Chi 


The requirements of axial symmetry and time independence imply 
Letra oi 4A h 
et Bog (x ,X ) where we have chosen x =¢ , x =T _ to be the angu- 
lar and time variables respectively. The gravitational field is invari- 
ant under the simultaneous transformation, $+ -d? , tT * -t since this 
leaves the angular velocity and thus the source of the gravitational 
field unchanged. This invariance of the metric requires 
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Greek indices are assumed to take the values 1 to 4, Latin indices 
Deon. 
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The coordinates on) can be chosen such that [1] 


B41 (ax*)? + 2g, dx" ae aa B59 (dx)? = eee es) aa (Ge) 
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Hence if we now define 
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“lee ee 
then equation (1.1.3) takes the form 
epee ae a) (ae de (ge edo 
- C2 (xt ,x7) (dt = v(xtyx*) dg)? (1.1.6) 


In order to obtain further simplification it is necessary to 
use the vacuum field equations. A variation of the argument used in the 


static case (see Synge [2]) leads to the result 
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where we have put x =p , x =z. It is convenient to define 
r 
C(p ,z) ac 3 
—) 
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where v and i) are functions of op and z. The final result is 
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where i , v,w are functions of op and z. 


The vacuum field equations [3] (i.e. Rug = 0 ) take the 
form 
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Equation (1.1.13) can be rewritten as follows: 
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The form of the above equation suggests the introduction of a function 


@ defined by 
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It has been demonstrated by Ernst [4] that the substitution 
Gaeta s pio Ginagite) 


implies that (1.1.12) and (1.1.13) are equivalent to the single complex 


equation 
Z 
(ReB)(VE) = (&) aes 
where V is the gradient operator for a three dimensional flat space. 


The imaginary part of (1.1.17) gives the following equation 


foro 
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which replaces equation (1.1.13). 
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To obtain the linearized equations for i , @,wv and y 
one assumes that all nonlinear terms, that is, terms involving products 
of field quantities or products of the derivatives of field quantities 


can be neglected. 


Thus we can write the linearized versions of equations (1.1.10), 
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The symmetries of physical systems often lead to the use of 
"preferred" coordinate systems, that is, coordinate systems in which the 
field equations take a particularly simple form. Since a disk is a 
degenerate oblate spheroid it is likely that oblate spheroidal coordi- 


nates will simplify the description of the disk. 


Oblate spheroidal coordinates are defined by [5] 
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where 0<e<™ , -l<n<i1 . These coordinates are ideally suited 


for a disk since the disk surface is given by e«=0. 
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The Laplacian of an arbitrary function F has the form 
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Separation of variables leads to the following solution to Laplaces 


equation 
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where P and On are respectively Legendre functions of the first 
and second kind. We shall later use expression (1.1.25) in the descrip- 


tion of the linearized field of the disk. 


it ane Linear zed Solution: 


The linearized equations written in terms of oblate spher- 


oidal coordinates take the form 
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To simplify the problem we shall make the following assump- 
tions: 
Ci) the angular velocity w as seen by an observer at 


infinity is constant. 
Cid) the gravitational and centrifugal forces on each parti- 


cle are equal. 


Since (ii) implies that the pressure is zero in the disk, we 


can write the solution to Laplace's equation (1.2.2) as 


A(x,y,Z) = - | el davdy | Ce oy) 
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where o(x',y') is the surface density of the disk and R is the distance 


from source point (x',y') to observation point (x,y,z) . 


Using expression (1.2.5) for A we can evaluate ale : 
the jump in the normal derivative of } due to the surface layer o . 


We obtain 
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Conditions (i) and (ii) may be expressed as 
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In terms of oblate spheroidal coordinates (1.2.10) has the 


form 
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As we observed earlier i takes the form 
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Condition (1.2.9) implies a= QO for m odd, and (1.2.11) implies 

aa = 0 for @y> 2. The coefficients ay and ay can be determined 
n 

from (1.2.8) and (1.2.11) and the requirement o be finite for (n=0) 


The result is 
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In the linear approximation the Bap are given by [2] 
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where Nye is the Minkowski metric and ae the energy momentum tensor 


of the matter distribution. 


According to the physical interpretation of the energy momentum 


tensor we have 
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From the above expressions we can determine the jump in the 


normal derivative of Bix ° 8 and Ben One concludes: 
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and in the linear approximation Be6 =w , so (1.2.18) leads to the 
result 
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From (1.2.3) one can now obtain the appropriate boundary 


condition for © , namely: 
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If we assume that © has the form (1.1.25) then the above conditions 


imply 
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At this point it seems convenient to introduce the dimension- 


less parameter 
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In this notation the Ernst function (1.1.16) can be written in 
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which can be regarded as the first few terms of a power series expansion. 
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It turns out that this point of view is useful in determining higher 


order approximations. 


To complete the specification of the linearized field we must 


determine ~ and v. According to (1.2.3) we have 
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which determines ~ up to an arbitrary constant. By the assumption 


of asymptotic flatness (i.e. € >» => p> 0 ) we obtain w in the form 
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where prime denotes differentiation with respect to the argument of the 


function. 


The solution to the v-equation (1.2.1) that is consistent with 


asymptotic flatness is 
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The metric (1.1.9) in terms of oblate spherical coordinates 


has the form 
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The linearized version of the above metric is 


2 
ds” = (1-24) [Aa?(e74n?) (42, + 


2 
1+ fie: 
: : (1.2.31) 


P(E Dah ea) tna edad eee CO) de evade 


More explicitly we conclude: 


Ore? 62 
g., 7 {1+ « 0° ies Se) Pod} 
I+e 
ee 
_ A (Ce +n_) NS mee 
on et prety (Lo-tyctome[QsePe =9Q57253) 
(42232) 
Byy 7A (te) (en) (1+ «0 1G, 2, - G PD 
sbalticecenactq) PO) 511 
bor ‘i oy roe is) ee 2 
Sepa 2 Z Pat Me nA'p! 
Er = Sor Te On ous (Cle ) (1-n )(QiPy ies 6 Q3P4)] 


where the terms inside the square brackets represent the 


deviation from flatness of the geometry. 
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CHAPTER II 


The Surface Energy Tensor 


I. Hypersurfaces in General Relativity [6]. 


Since we are concerned with a disk of infinitesimal thickness 
in a four-dimensional space-time we shall give a short discussion of 


the geometry of hypersurfaces. 


In order to completely describe a hypersurface it is necessary 
to specify the intrinsic properties of the hypersurface (i.e., those 
properties which are independent of the imbedding space), and the so- 
called extrinsic properties which describe how the hypersurface is imbed- 


ded in the surrounding space. 


The extrinsic geometry of a hypersurface can be described by 


a 3x3 symmetric tensor which is defined as follows. 


Suppose we are given a four-dimensional Riemannian space-time 
V , in which is imbedded a smooth hypersurface ) . Further suppose 


that ) is defined by 


neue (en) ies (2.1.1) 
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The ¢ can be considered as intrinsic coordinates on the 


hypersurface. 
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An infinitesimal displacement in )} 


dx* =f det ORs 
i 
d€ 
defines the quantity f° which transforms as a contravariant four- 


ai 


vector for fixed i , and as a covariant surface-vector when a is 
’ ; a 
fixed. Also the triad of four-vectors (f £ 


linearly independent and tangent to i) 5 


ie n> is a space-like unit vector normal to ) then we note 


that 
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is (for each fixed i ) a four-vector tangent to the surface. Hence 


there exists quantities Ke such that 


‘cs (2.1.5) 


I 

Aw 
as 

Fh 


0) 
on 
6€., 

i 

K.. is called the extrinsic curvature tensor of the hypersurface. 
1j 


Since the surface metric tensor is given by 
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The extrinsic curvature tensor can be used to define a 


singular hypersurface or surface layer in general relativity as follows. 
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For a justification of the name surface energy tensor we 


refer the reader to Israel [6]. 


The tensor S44 defines three linearly independent eigen- 
vectors one of which is timelike. Hence the surface density o can 


be defined by the eigenvalue equation 
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II. The Surface Energy Tensor of the Disk. 


We now apply the previous results to the special case in 
which ) is the surface of a disk. In oblate spheroidal coordinates 
the surface is defined by e¢ = 0. We will use the following convention 


for the coordinates 
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(Subscripts on functions indicate partial differentiation with respect 


to that subscript.) 
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The above explicit expressions for the surface energy tensor 
for the linearized gravitational field can be used to determine the 


surface rest mass density. According to equation (2.1.10) we must solve 
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CHAPTER IIT 


Higher Order Approximations 


I. Outline of Method. 


In chapter I we noted that the Ernst function appeared to 
have a natural expansion in terms of the parameter a. By choosing a 
particular form of the Ernst function it turns out that with very little 


effort the second order approximation can be obtained. 


Recall that in chapter I, equation (1.1.16), we introduced 


the Ernst function 


It is advantageous to rewrite the above expression in terms of two real 


valued functions A(e,n) , Ble,n) as follows 
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These equations are ideal for approximation purposes since A 


and B are harmonic to orders ie and ae respectively. 
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Recall that w is determined by the equation (1.1.15) which 


in oblate spheroidal coordinates take the form 
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Finally in order to complete the description of the gravitational 


field it is necessary to calculate wv. 
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II. The Gravitational Field to Order ae 


As we noted earlier knowledge of the linearized gravitational 
fieid, ise. A, ; BY is sufficient to determine the gravitational field to 
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where the constant of integration is zero by the requirement of asymp- 


totic flatness. 


The functions oF and Vo which correspond to i =a by 
and = ay, are obtained by the integration of (1.1.23) and(_1.1.24). 


The result is given in appendix A. 
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#t1L. The Surface Energy Tensor to Order o - 


The functions  , w and wv can now be used to calculate the 


surface energy tensor by the technique used in chapter II. 


Exact expressions for the surface energy tensor of the disk 


are given in equation (2.2.8). Neglecting all terms of higher order than 
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0, one arrives at the following expressions: 
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We note that a = 0 when n=0 i.e. on the edge of the 
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disk. Also the angular momentum density S goes to zero at the 
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center (n=1) 


IV. Physical Properties of the Disk. 


The metric and the energy tensor can be used to calculate 


several physical parameters which characterize the disk. 


The gravitational mass m of the disk can be calculated from 


the following integral [7] 
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From figure 1 we observe that the surface mass density decreases 
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from its linearized value with the maximum decrease at approximately jn = 
Bardeen [8], in his discussion of a uniformly rotating disk, has performed 
numerical calculations which show that for large expansion parameters, the 


proper surface density has a maximum near the edge of the disk. 
The angular velocity (2.2.15) as seen by an observer at infinity 
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It is evident that w attains its maximum at the rim of the 


disk and its minimum at the center. 


The surface energy density 06 can be used to obtain the rest 


mass [9] 
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Making use of (3.4.2), (3.4.7), (3.4.10), (3.4.8) and (3.2.4) 


the rest mass (3.4.9) is 
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is defined as the energy required to remove matter from the disk to 
infinity, i.e. the energy required to overcome gravitational forces 


and rotation, thus 
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Thus energy is released when the disk is formed. 


Using (3.4.11) and (3.4.3) we can determine the fractional 
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According to Bardeen's numerical computations the fractional 


binding energy is a positive monotonically increasing function of the 


expansion parameter. 
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CHAPTER IV 


Applications of the Metric 


I. Geometry of the Disk. 


Since the gravitational potentials 8,4 are assumed to be 
a] 


continuous functions, it follows that the exterior field uniquely 


determines the metric on the disk. 


In special relativity an element, df of spatial distance is 
defined to be the interval between two infinitismally separated events 
occurringat the same time. However in general relativity proper time 
is a function of position. Hence, it is necessary to make a correction 


to the t = constant form (see Landau and Lifshitz [10]) 
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that of flat space in oblate spheroidal coordinates for n approaching 


+ 1 where the radius of the disk is given by 


' = — 
A A exp ( ay aan? ; C43) 
Thus we rewrite (4.1.2) as 
i eee 
ae? = CAB eee eyhia2, oe 
Sh 
=X 2(A_-A) 
SoCo 1 Gees. Se ecte (ies) 


where Ms is the i - function evaluated at the center of the disk 
(n = +1) . The bracketed quantities above give the deviation of the 
disks geometry from the flat space determined by elementary flatness, 
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In the above expression rotational effects are observed in the 
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symmetry for constant n is 
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II. Asymptotic Properties. 


We would like to examine the disk metric for large e¢« . Since 
we have postulated asymptotic flatness it is necessary that ¢ = 
reduces the disk metric to that of Minkowski space-time. If we take 
into account first order corrections from flatness (ieneei) we expect 
a certain simularity to the asymptotic form of the Schwartzschild metric. 
If we consider rotational terms Cie eere) the far field of the disk 


should have the form of the linearized metric of a rotating ball of fluid. 


In order to substantiate the above expected properties we will 
compare the disk metric for large e to the asymptotic form of the Kerr 
metric [11] which is the only known stationary exact solution which obeys 


asymptotic flatness, 


The Kerr metric can be derived by the Ernst function [4] tech- 
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and we have set A=1 in (4.2.1). 


We note that the relationship between polar coordinates, oblate 


spheroidal coordinates and prolate spheroidal coordinates is given by 
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Expanding the above functions for large x and rewriting the 


expansions in terms of (r,6) gives the following results: 
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Comparing the Kerr metric (4.2.5) to the disk metric 4.5226) 


we observe that 
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must be the graviational mass of the disk, which agrees with( 3.4.3). 


The parameter a of the Kerr metric must correspond to 
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This is the value of a which could be obtained by inter- 
preting the quantity (-ma) of the Kerr metric as the total angular 


momentum of the disk given by (3.4.5). 


Rewriting the disk metric in terms of the above two parameters 
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The Schwartzchild metric can be obtained from the Kerr metric 


by letting the parameter a go to zero. Comparing (4.2.9) atid (4.2.5)) 


wa 
» . 7 
> 
tc 
’ 
i 
* 
{ 
4 
eu 
- 
‘4 


aaiP 


SiWs5m. 1704 


‘ 
(88.508) bbs ad gatas sce eee oF 
- 7 : 5 > 


> 7 
7 
i 
a a . 4 
ae ‘ . ' v 4 e. - my *» 
jou Wakb ett oF (3) S) oe ee) Geet, 
: - 5 
er ee 
A - , *ec3 | oe ys 
; 
4 on 
i ‘ ie : i 
a in 7 a 
‘J sly 
i * : 
| j a ‘Beer 
i 
Ey yx ‘ 4 — mein wm b, 


t f - 
r vv b 
. 
i. o, 
»% 
, 4 ‘ Dy my 
? 
A 
P Pea, 
| a bs ro ey : ry i fe i f Seth 
! 7 - 
7 
e 
ay rail 
j ie 
thea 
ow Orr n 
‘ re . 
~— tu f - —. ism = oe 
| a f #* ss 
~~ he 
of \ 
AS h 
» 
iY 
a : 
é r 
. wt a = 
* ie a z 
> be { ' . t wid - | Seeded = 
ope a 
ge ies 
i 
€ oc ; 
- i n= : i > 
SE cot! Le Te eee tae = & 
f ee ‘oc amel | ee ad ad a es 1 ee, 
ey r ’ 


i] 
, % i 


uJ ih 


7 ; = | . 
aris mot? bagkeida sd nen oabris 


j a 


aa 


e 
7 a 


pale adi 
: iow 


BE 


with a= 0 shows the disk deviates from the Schwartzchild metric in 


1 P fom 
the (=) terms with the minimum deviation obtained for an observer 
r 


along the axis of symmetry. If we consider the effects of rotation 

(a#0) the two metrics coincide to order SS) . Since the (=) terms 
of the Kerr metric in Schwartzchild-like eben Sess are the same as the 
linearized metric of a rotating ball of fluid [12] we infer that the disk 


metric is asymptotic to the linearized metric of a rotating ball of fluid. 


For a further discussion of the disk metric and the Kerr metric 


in the strong field limit we refer the reader to Bardeen [8]. 


III. The Gravitational Red Shift. 


In this section we obtain an expression for the change in the 


frequency of a photon due to the gravitational field of the disk. 


Consider a radiating atom at the center of the disk (0,1) 
and an observer at (e,n) - The proper time intervals at the center and 
at (e,n) are denoted respectively by 
i 
DN os CET te: 


(4.3.1) 


If an atom at the center (in proper time On: emits n-waves 
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At the observation point the receiver would observe the same n-waves in 
some proper time interval At . Evidently the observed frequency vy 


and the proper time At would be related by 


n= v At e (47333) 
Hence, 

AT bet 5 
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The red shift is given by 
fet S 
" Set 


From (3.2.4) we have 
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Using (4.3.6), (4.3.5) and the appropriate expressions in appendix B we 
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obtain (correct to order a ) 
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Evidently the maximum red shift occurs for an observer at 


infinity where r, = x =O). Im this»casen(4.3.7)) reduces: to 
Av te Teo ee ae Oe 27k 1, ees 
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On the disk (e=0) the red shift (now a function of 7 ) is 


given by 
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From the above expression it is observed that to order a 
the red shift is identical to that of a static disk. The rotational 


term in the a approximation decreases the frequency of light emitted 


at the center of the disk. 


The maximum red shift on the disk is obtained at the edge. 
Setting n= 0 we obtain 


Z PLY ous 2 
a 3k + at oan 6 Zia 25k CeO) 


y sec2s 512 a fe 24,576 36 ] 


o | n=0 
e=0 


To obtain the red shift of a photon emitted at the edge of the 


disk and received at infinity we subtract (4.3.10) from (4.3.8) which gives 
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If we only consider fe - terms, half of the total energy loss of the 
photon in the z= 0 plane is obtained when the photon reaches the 

edge of the disk the other half is lost in leaving the external gravita- 
tional field of the disk. We also observe that rotation doesn't effect 
the red shift of a photon emitted at the rim of the disk in the order of 


approximation we have considered, 


The red shift of a photon emitted at the surface of a spherical 


mass distribution of mass m and radius A is given by 


et (4.3.12) 


2 3 
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it 3 at ta @® 3,072 ‘A (4.3.13) 


Thus the energy loss will be slightly greater for a photon leaving a 


disk in the z= 0 plane than a spherical object. 


IV. Lense-Thirring Effect. 


In this section we consider the effect of the component 


beg 
of the metric tensor which gives rise to the so-called “dragging of the 


inertial frames'' of an observer. 


Hartle and Sharp [17] have shown that the axially symmetric 


stationary metric can be written in the form 
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da = — y“dt* 4 Q*dr* + r°x7 fae? + sin?e (dg - L dt)2] (resi) 


where H,Q, K and L depend on (r,6). By an argument analogous to that 


given by Hartle [17] the function L is the angular velocity of an observer who 
falls from infinity to the point (r,6) . Thus L is the rate of rota- 


tion of the inertial frame with respect to flat space, 


If (1.1.9) and (4.5.1) are compared it follows that 
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A? (+07) (1-n") 
On the surface of the disk (4.5.3) becomes 


3 5k 1 (143n°) 


° 4.5.4 
32A ( ) 


L = a 


We observe the rotation is greatest at the center of the disk. 


Now the w of (3.4.8) is the angular velocity of the fluid as 


seen by an observer at infinity. 


Hence, 
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is the rate of rotation of the fluid with respect to a local inertial 


frame as seen by an observer at infinity. 


From (4.5.5), (4.5.4) and (3.4.8) one concludes that for points 


on the disk 


BOO 0 on 


oy Z 

W 6A 06A (3+25n ) e (4.5.6) 
Thus the rate of rotation of the fluid with respect to an inertial frame 
is less than with respect to an observer at infinity. The decrease is 


observed to be maximum at the center and minimum at the edge, 
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CHAPTER V 


Summary and Conclusions 
ee) eee One rus LOns) 


In this final chapter we would like to summarize the results 


that have been obtained. 


The first chapter was concerned with determining the linear- 
ized field of a uniformly rotating, infinitismally thin disk. The choice 
of uniform rotation was governed by two considerations (i) the work of 
Hunter [13] (a Newtonian treatment of infinitismally thin disks) has 
shown that the simplest possible mass distribution o = 2 oon leads to 
uniform rotation, (ii) the possible use of a thin uniformly rotating 
disk as a model for a galaxy (rotation is observed to be approximately 


uniform in certain galactic systems [14]). 


Concerning the astrophysical implications of a thin rotating 
disk it should also be emphasised that Hunter has also shown that an 
infinitismally thin Newtonian disk is unstable. Since general relativity 


enhances instabilities we would not expect the relativistic disk to be 


stable. 


There are two possible generalizations which may improve the 
stability of the disk, namely, (i) differential rotation, (ii) assume 
the disk has thickness. Differential rotation can be obtained for a disk 
by assuming o is a polynomial in n . The techniques used to determine 


the linearized metric could then be carried out. However, the work 
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involved would be considerable. For a further discussion of stability 


we refer the reader to Bardeen [8]. 


The second chapter contains explicit expressions for the energy 
tensor. It is shown that the proper surface density and angular velocity 
(correct to order 0°) correspond to their appropriate Newtonian values. 
From equation (2.2.12) it follows that rotation does not effect the Newton- 


? : : 4 
ian surface density until one considers terms of order a . 


In the third chapter the gravitational field was determined to 
order ee - Integral expressions for the functions which determine the 
8 
gravitational field were obtained correct to order a . However, the 


integrals are so unwieldy that these results are not included. 


It was also noted that the angular velocity w , (3.4.8) is no 
longer constant but becomes a function of mn. The matter rotates faster 
at the edge than at the center, which is opposite to what one would want 
for a realistic galaxy. In order to obtain a disk which is uniformly rota- 


ting one would have to choose a different initial mass distribution. 


We also observed that the calculated binding energy and fractional 
binding energy are positive which is in agreement with that of a incompres- 


sible fluid sphere. [15]. 


Chapter IV contains a discussion of certain properties of the disk 
metric. In particular we considered such effects as the gravitational red 


shift, asymptotic properties, and the dragging of inertial frames. 
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When the metric of the disk was compared to the Kerr metric 


the observed deviations were in the terms of order 4) . The devia- 
rs 


tion was minimal for an observer on the axis of symmetry. 


The red shift which we considered was due only to the gravita- 
tional field. This is not the red shift which is physically observed since 
there is an additional shift in the frequency of light due to the rotation 
of the disk. The effects of the rotational graviational field on the 
photons energy loss appeared only in the a approximation causing a 


greater red shift than in the static case. 


The effects of rotation were also observed in the so-called 
, 3 
"dragging of the inertial frames". This effect was calculated to order a 
It was concluded than an inertial observer sees the fluid rotate faster at 


the edge than at the center. 


By determining the gravitational field of a rotating disk we have 
also obtained the field to a static disk. To see this we note from equations 
(1.1.10), (1.1.11) and (1.1.12) that »=0 gives equations(1.1.20) which 
we have solved. Thus, although we can't let be zero in our above solution 


the results obtained for = 0 is the same as the results which one would 


obtain by solving the static field equations. 


In conclusion we would like to comment on the relative strengths 


: 2 271m 
of the above effects. The expansion parameter is given by a = TSA 


If one considers the galaxy NGC 2782 which obeys a linear rotation law [14] 


and has a mass and radius given by 
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m= 17 x 107° ma R= 1.3 kip.e¢ 


it follows that 


We thus observe that for a normal galaxy our expansion para- 
meter is indeed less than one. One might therefore expect that the 
description of such galaxys in terms of the linearized solution to be 
quite good. Of course for the more interesting case of strong fields 


(i.e. where M ~A) the above analysis would be inadequate. 


Finally we refer the reader to Cameron [16] for a detailed 
discussion of the astrophysical applications of rotating disks and to 


Bardeen [8] for a complete numerical treatment of the disk problem 


including the strong field limit. 
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APPENDIX A 


Integration of (1.1.23) and (1.1.24) using do and Vo 


given by (1.2.13) and (1.2.28) gives the following result for 
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The above expression satisfies both asymptotic flatness and 


elementary flatness (v>0O as n>*tHt1). 
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APPENDIX B 


We list the functions A ; hy A i - Vy >» Vz» Vo and their 


derivatives with respect to e evaluated on the disk (e=0). 
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